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The effect of craters on the lunar neutron flux 

V. R. Eke 1 , K. E. Bower 1 , S. Diserens 1 , M. Ryder 1 , P.E.L. Yeomans 1 , L. F. A. 
Teodoro 2 , R. C. Elphic 3 , W. C. Feldman 4 , B. Hermalyn 5 , C. M. Lavelle 6 , D. J. 
Lawrence 6 


The variation of remotely sensed neutron count rates is 
measured as a function of cratercentric distance using data 
from the Lunar Prospector Neutron Spectrometer. The 
count rate, stacked over many craters, peaks over the crater 
centre, has a minimum near the crater rim and at larger 
distances it increases to a mean value that is up to 1% lower 
than the mean count rate observed over the crater. A simple 
model is presented, based upon an analytical topographical 
profile for the stacked craters fitted to data from the Lunar 
Orbiter Laser Altimeter (LOLA). The effect of topography 
coupled with neutron beaming from the surface largely re¬ 
produces the observed count rate profiles. However, a model 
that better fits the observations can be found by including 
the additional freedom to increase the neutron emissivity of 
the crater area by ~ 0.35% relative to the unperturbed sur¬ 
face. It is unclear what might give rise to this effect, but 
it may relate to additional surface roughness in the vicini¬ 
ties of craters. The amplitude of the crater-related signal 
in the neutron count rate is small, but not too small to de¬ 
mand consideration when inferring water-equivalent hydro¬ 
gen (WEH) weight percentages in polar permanently shaded 
regions (PSRs). If the small crater-wide count rate excess is 
concentrated into a much smaller PSR, then it can lead to 
a large bias in the inferred WEH weight percentage. For in¬ 
stance, it may increase the inferred WEH for Cabeus crater 
at the Moon’s South Pole from ~ 1% to ~ 4%. 


1. Introduction 

Cosmic ray interactions with planetary surfaces lead to 
nuclear fragments being released in the regolith. The study 
of neutrons that avoid nuclear recapture and subsequently 
escape through the surface provides a route to determining 
the abundance of various nuclei near the surface of those 
bodies [Lingenfelter et al, 1961; Metzger and Drake, 1990; 
Feldman et al., 1991]. Of particular interest is the epither¬ 
mal neutron flux (energies in the range 0.3 eV< E < 0.5 
MeV), because of its sensitive dependence on the hydrogen 
abundance in the top ~ 70 cm of regolith [Feldman et al., 
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2000]. The first experiment to search for lunar hydrogen in 
this way was the Lunar Prospector Neutron Spectrometer 
[LPNS, Feldman et al, 2004]. Feldman et al. [1998] found 
that there were polar dips in the epithermal neutron count 
rate, implying the existence of polar near-surface hydrogen. 
Furthermore, the lack of a corresponding feature in the fast 
neutrons with energies exceeding 0.5 MeV [Feldman et al., 
1998] or the thermal neutrons with E < 0.3 eV [Lawrence 
et al., 2006] suggested that any hydrogen-rich layer of ma¬ 
terial should be buried beneath 5—10 cm of hydrogen-poor 
material. 

The omni-directional LPNS, when orbiting at 30 km, had 
a spatial footprint with a full-width half-maximum (FWHM) 
on the lunar surface of 45 km [Maurice et al., 2004]. In or¬ 
der to suppress statistical noise, Feldman et al. [1998] and 
Feldman et al. [2001] binned the data into ~ 60 km x60 
km pixels. This is large relative to the sizes of most perma¬ 
nently shaded regions. However, Eke et al. [2009] showed, by 
stacking data and using a pixon-based image reconstruction 
technique to improve the spatial resolution while suppress¬ 
ing the noise, that these count rate dips could, in a statistical 
sense, be associated with the permanently shaded regions. 
This result was confirmed by Teodoro et al. [2010] using a 
more accurate set of permanently shaded regions defined 
by the SELENE laser altimeter [Noda et al, 2008]. The 
count rate dip inferred for Cabeus crater corresponded to a 
(1 ± 0.3)wt% water equivalent hydrogen (WEH) according 
to the regolith composition model of Lawrence et al. [2006], 
which has a semi-infinite layer of ferroan anorthosite (FAN)- 
type soil with varying amounts of H 2 O. 

Another experiment, the Lunar Exploration Neutron De¬ 
tector (LEND), contained sensors called the Collimated 
Sensor for EpiThermal Neutrons (CSETN) and the Sensor 
for EpiThermal Neutrons [SETN, Mitrofanov et al., 2010]. 
These mapped the Moon from the Lunar Reconnaissance 
Orbiter at an altitude of 50 km, ~ 20 km above the or¬ 
bit of Lunar Prospector. Thus, one should not expect the 
SETN instrument to provide competitive results relative to 
the LPNS. Furthermore, comprehensive analyses of the data 
returned from the CSETN have demonstrated that the col¬ 
limator did not perform well enough to fulfil its mission ob¬ 
jectives [Lawrence et al, 2011a; Miller et al, 2012], with 
the vast majority of lunar neutrons being uncollimated [Eke 
et al., 2012] and an effective FWHM much larger than that 
of the omni-directional LPNS [Teodoro et al, 2014]. In view 
of the difficulties associated with the interpretation of this 
data set, these data will be considered only briefly in this 
paper. 

The Lunar Crater Observation and Sensing Satellite 
(LCROSS) impacted into Cabeus crater in 2009 and the 
resulting ejecta plume was analysed to give a value of 
(5.6 ± 2.9)wt% WEH [Colaprete et al., 2010]. While sta¬ 
tistically consistent with the LPNS result, the most prob¬ 
able value is over five times the LPNS-inferred value. The 
reanalysis of the LCROSS data by Strycker et al. [2013], 
which gave (6.3 ± 1.6)wt% WEH, is inconsistent with the 
LPNS result. These comparisons would be affected if the 
hydrogen detected by the LPNS was not uniformly spread 
across the surface within the large resolution element, which 
is approximately 1000 times as long as the crater produced 
by the LCROSS impact [Schultz et al., 2010]. The LPNS 
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and LCROSS results sample somewhat different depths into 
the regolith. Thus, any variation in hydrogen content with 
depth could also lead to a difference between the hydrogen 
abundances inferred from the two separate methods. One 
assumption that is implicit in the studies of craters using 
the LPNS data is that neutron count rates are not explicitly 
affected by the surface topography. A new model will be 
presented in this paper to quantify the topographical effect 
from craters on the neutron count rate. 

The Chandrayaan-1 M 3 results interpreted as implying a 
particular excess of water or hydroxyl molecules in Gold¬ 
schmidt crater [ Pieters et al, 2009] prompted Lawrence 
et al. [2011b] to re-examine LPNS data in this region in the 
context of a two-layer regolith model, with the surface layer 
being hydrogen-rich. This contrasted with previous Monte 
Carlo modelling of the lunar regolith, where the hydrogen 
had been buried under a dry layer of regolith [ Lawrence 
et al., 2006]. After removal of the trends caused by bulk 
composition, the thermal and epithermal data in the vicin¬ 
ity of Goldschmidt crater were compared with the models 
to investigate the sensitivity of neutron measurements to 
the depth distribution of hydrogen. Lawrence et al. [2011b] 
concluded that it was necessary to understand more about 
systematic variations at the 1 — 3% level before definitive 
conclusions could be reached. If crater topography does pro¬ 
vide small systematic variations in neutron count rates then 
it needs to be understood in order to progress. 

When studying the LPNS count rate, Feldman et al. 
[2001] found that ‘local maxima overlay the floors of large, 
flat-bottomed craters’. They did not speculate as to what 
this implied, but the possibility of a topographical effect on 
the measured neutron count rate is one that could create a 
systematic bias in the values of WEH inferred above per¬ 
manently shaded polar craters. To date there has been no 
systematic, quantitative study of the imprint of topographi¬ 
cal features on the detected orbital neutron count rate. It is 
important to quantify the impact of topography on the emit¬ 
ted lunar neutron flux because many of the results from the 
LPNS involve small changes in count rates measured over 
craters. 

Section 2 contains a description of the neutron and to¬ 
pography data being used. Fits to the crater average topog¬ 
raphy are given in Section 3. The variation of neutron count 
rate as a function of distance to the crater centre is shown 
in Section 4, for a variety of different subsets of craters. In 
Section 5, a simple model is presented for how the neutron 
count rate changes as a function of detector distance from 
the crater centre. This model is confronted with the data, 
and the implications for our understanding of the regolith 
are investigated. Section 6 discusses the implications of this 
work for quantitative estimates of cold-trapped hydrogen, 
and conclusions are drawn in Section 7. 

2. Data 

Maps of the lunar neutron count rate, a set of predeter¬ 
mined lunar craters and a digital elevation map are neces¬ 
sary to calculate the neutron count rate profiles near craters. 
The data sets to be used here, which are all available from 
the Geosciences Node of NASA’s Planetary Data System 
(PDS 1 ), are described in this section. 

2.1. Lunar Prospector neutron data 

The Lunar Prospector spacecraft spent one year at 100 
km altitude, then seven months at 40 — 30 km. PDS time 
series data from the thermal, epithermal and fast neutron 
detectors, processed as described by Maurice et al. [2004], 
are used in this study, with the focus mainly on the low- 
altitude subset. Some results from the high-altitude period 
will also be shown for comparison, but the default choice is 


to consider only data for which Lunar Prospector was at an 
altitude less than 45 km. Using different energy neutrons 
is desirable because of their differing responses to changes 
in regolith composition. Also, the thermal neutrons probe 
further into the regolith than the epithermals, whereas the 
fast neutrons typically sample nearer to the surface than the 
epithermals. 

2.2. Lunar Exploration Neutron Detector data 

Data from the first 15 months of the mapping orbits are 
used for both the LEND SETN and CSETN detectors, to 
compare with the results from the LPNS. For the CSETN 
measurements the background due to cosmic rays striking 
the spacecraft itself is removed statistically following the 
procedure described by Eke et al. [2012], The remaining 
count rate is comprised of two distinct lunar components, 
where the detected neutrons originate either from within or 
outside the collimator’s geometrical field-of-view. One can¬ 
not determine from which component individual neutrons 
originate. 

2.3. Crater list and topographical data 

The list of craters produced by Head et al. [2010] from the 
Lunar Orbiter Laser Altimeter (LOLA) topographical data 
is used. This consists of 5185 craters with radii of at least 10 
km, distributed over the entire lunar surface. In this study, 
various different selections of craters are made, based on the 
radii, r c , and the central locations given in this list. The 
variable r will be used here to represent arclengths along an 
unperturbed spherical surface, whereas the variable x rep¬ 
resents the distance from the symmetry axis (z) of a crater. 
Thus, the measured crater diameters are really 2x c in this 
nomenclature, x and r are related via 

x = r m sin ( — ) , (1) 

VW 

where r m = 1737.4 km is the lunar radius. This equation 
implies that x c will be within 0.1% of r c for crater radii less 
than 100 km, so the variables x c and r c will be assumed to 
be equal for the rest of this study. The angle subtended at 
the lunar centre by the crater radius is 



The global topographic map from LOLA [Smith et al., 
2010] with (1/64)° resolution is used to measure the crater 
topographical profiles. This corresponds to ~ 0.5 km reso¬ 
lution at the equator, which is more than sufficient for the 
approximate modelling of crater topography as a function 
of crater radius that is necessary for the neutron count rate 
model presented in Section 5. 

The epithermal neutron count rate measured by the om¬ 
nidirectional LP detector changes by approximately 10% 
across the whole Moon. This variation is dominated by 
known changes in regolith composition. Any systematic to¬ 
pographical effects are expected to be at the level of ~ 1%, 
as noted by Feldman et al. [2001]. This anticipated varia¬ 
tion is sufficiently small that it is necessary to stack together 
craters of similar size in order to reduce the statistical uncer¬ 
tainties. In addition, the stacking averages away azimuthal 
anisotropy that exists in the crater sample, making radial 
profiles an appropriate way to represent the results. To pro¬ 
duce a more homogeneous set of craters to stack together, 
both in terms of topography and composition, only craters 
in ‘highland’ regions will be considered in this study. This 
means only craters on the far side of the Moon and with 
latitudes greater than —20° will be included in the stacking 
procedure. These cuts leave just 2216 craters. This choice is 
important for some of the results presented later involving 
thermal and fast neutrons, which are both more sensitive 
than epithermal neutron fluxes to iron and titanium abun¬ 
dances. 



EKE ET AL.: LUNAR CRATERS AND NEUTRON FLUXES 


X - 3 



Figure 1. Model crater profile (bold line) that is fitted 
to the LOLA topographical data for the selected craters. 
The model takes into account a flat infill region at radii 
less than Xi , a spherical cap depression out to x c and an 
outer slope, uplifted at the crater rim by an amount u, 
that returns to the unperturbed surface at x e . The cen¬ 
tre of the Moon is used as the origin of the coordinate 
system. 


the gradient of the outer slope is given by 

r m cos 9 e - (r m cos 6 C + u) . . 

9= ---• (6) 

X e — X c 

The height of the crater relative to an unperturbed sur¬ 
face can be found using 

h(x) = \Jx 2 + z(x) 2 - r m , (7) 

where 2 is defined as zero at the lunar centre. For x > x e , 
h = 0 for the model crater. At x = 0, 

z(x = 0) = 2o = r m cos# c — d T u + di. (8) 

At radii where the presence of the crater perturbs the 
surface, the height can be inferred using the following ex¬ 
pression for 2 : 

! 2 o if x < Xi ; 

2 0 - di T r curv - Vrc UIV - x' 2 if Xi < x < x c ; 

2 C -I- g(x — x c ) if x c < X < x e ■ 

f _ ( 9 ) 

2 C represents z(x c ) = zo — di + r cur v — \A1Lrv — 

A least-squares minimisation was performed to find the 
best-fitting sets of the parameters [d,Xi,u,x e ] for the sub¬ 
sets of craters of different radius. For all subsets of craters 
with radii of at least 20 km, the region at x < 0.2x c was 
excluded from the fit, because a central peak, not included 
in the model, often exists. The midpoint of the crater radius 
bin was chosen for the x c of the model to be fitted to the 
stacked profile. 


3. Crater topography 

The model for the topographical effect on the neutron 
count rate described in Section 5 needs to assume a particu¬ 
lar crater profile. This section first describes the functional 
form of the assumed crater profile, then measures it using 
the LOLA digital elevation map by stacking radial profiles 
for craters in the chosen range of sizes. 


3.1. Model crater profile 

Rather than having a general topography, the model 
craters in Section 5 are considered to have azimuthally sym¬ 
metric profiles of a kind shown in Figure 1. These consist of 
a spherical cap depression of depth d measured down from 
the x — y plane containing the crater rim, with a central, 
flat ( dz/dx = 0) infill region extending out to a radius Xi, 
and with a maximum depth, at x = 0 of di, where the i 
subscript refers to the infill region. The radius of curvature 
for the spherical cap part of the crater is then 


_ Xc + d 2 

~ 2d 


( 3 ) 


and the maximum infill depth, measured from the base of 
the spherical cap to the infill surface, is given by 


di =r 


curv 



( 4 ) 


This crater is uplifted parallel to the crater axis (the 2 di¬ 
rection) by a distance u, with an outer slope of constant 
gradient, dz/dx = g, back to the unperturbed surface at a 
perpendicular distance of x e from the symmetry axis of the 
crater. Defining 


8 e = sin 


Xe 

r m 


( 5 ) 


3.2. Crater topography fits 

For each crater, digital elevation map measurements 
within 3r c were used to construct the relative height profile 
as a function of r/r c , where r represents the arc length from 
the crater centre to the spacecraft nadir. The zero of height 



r/r c 


Figure 2. Mean radial LOLA topographical profiles for 
craters with radii in the ranges 10 — 20 km (black filled 
circles) and 40 — 50 km (blue open circles). Statistical 
errors on the mean profiles are smaller than the symbol 
sizes. Solid lines show the least-squares model fits to the 
data sets as detailed in equations (1-9). 
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Table 1. Least-squares parameter values for the crater topography model fitted to the stacked LOLA data as a 
function of crater radius. 3. 


Crater radius, r c /km 

Number of craters 

Depth, d/km 

Infill radius, x\/x c 

Edge of outer slope, x e /x c 

Uplift, u/ km 

10 - 20 

1264 

2.00 

0.30 

2.1 

0.45 

20 - 30 

482 

2.50 

0.42 

2.1 

0.60 

30-40 

219 

2.75 

0.48 

2.0 

0.70 

40 - 50 

115 

3.25 

0.56 

2.2 

0.85 

50 - 60 

47 

3.70 

0.66 

2.1 

1.15 

60 - 80 

42 

3.40 

0.70 

1.8 

1.00 

40 — 50 deep 

57 

4.05 

0.48 

2.2 

1.15 

40 — 50 shallow 

58 

3.05 

0.69 

2.1 

0.65 


3 The depth, fractional infill radius, fractional edge of outer slope and uplift are found on grids with resolution 0.05 km, 0.01, 0.01, 
and 0.05 km respectively, which are larger than the statistical uncertainties on these parameters. The midpoint of the crater radius 
range is used to calculate the parameters for each stacked profile. 



12 12 12 


r/r c r/r c r/r c 


Figure 3. Stacked, normalised radial LP neutron mean count rate profiles for different radius craters. 
Thermal, epithermal and fast neutron results are shown in the left, centre and right panels respectively. 
The different colours correspond to craters in the radius ranges 10 — 20 km (red), 20 — 30 km (magenta), 
30 — 40 km (orange), 40 — 50 km (black points), 50 — 60 km (green) and 60 — 80 km (blue). A black dotted 
line represents the epithermal neutron profile for 40 — 50 km craters when the time series is sampled in 
32 s observations, like the fast neutron data set, rather than 8 s like the thermal and epithermal time 
series. 


for each crater is defined as the mean height in the range 
2.5 < r/r c < 3. Each measurement provides an estimate of 
the relative height at its r/r c - The statistical uncertainty 
on the estimated mean height is just the square root of the 
ratio of the variance of the individual measurements within 
a given bin in r/r c to the number of observations in that 
bin. The craters were stacked by crater radius, because the 
typical crater shape varies systematically with crater radius. 

The crater set with 40 < r c /km< 50 was further subdi¬ 
vided by depth to see how this affected the neutron count 
rate profiles. To split the crater subset by crater depth, in 
order to investigate the effect on the neutron profile, the 
depth of each crater is defined as the difference between 
the average heights in the radial ranges (0.95 — l)r c and 
(0.2 — 0.3)r c . The central region is once again avoided to re¬ 
duce any systematic effect due to central peaks. While this 
statistic might, in some instances, reflect subcraters rather 
than the larger scale topography, it at least serves as a sim¬ 
ple way to separate deep and shallow craters with the same 
radius. 

Figure 2 shows the azimuthally averaged mean radial to¬ 
pographical profiles for craters of different sizes, as measured 
using craters from the Head et al. [2010] list and LOLA to¬ 
pographical data. Statistical errors on the mean profiles are 
smaller than the symbol sizes. It is apparent that for the 
larger craters there are central peaks that are not included 
in the model profile, as described in the previous section. 


The best-fitting models are also shown in Figure 2, from 
which it can be seen that the model becomes increasingly 
inappropriate for larger craters. A flat, dz/dx = 0 central 
region does not translate to a constant height relative to 
the unperturbed spherical surface, which would provide a 
better fit to the r c > 70 km craters. Also, the constant 
dz/dx outer slope, at large distances, can lead to h < 0 on 
the outer uplifted slope; a feature not present in the obser¬ 
vations. Despite these shortcomings in the model, it does 
capture the main features present in the measured average 
topographical profiles, and the extent to which the model is 
inadequate is not quantitatively significant for the neutron 
count rate results in subsequent sections. 

Table 1 lists the best-fitting model parameters for a set 
of different crater size ranges. These values are used in Sec¬ 
tion 5 for the model predicting the topographical effect on 
the neutron count rate profiles observed by the orbiting de¬ 
tector. The depth parameter, d, only represents the depth 
of the crater when there is no infill so, as can be seen by 
comparing the values in Table 1 with the data in Figure 2, 
the actual crater depths from rim to minimum are typically 
much smaller than d. This is particularly true for the larger 
craters, where the best-fitting infill region extends to a larger 
fraction of the crater radius. 

4. Neutron count rate profiles 
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For each crater, time series observations within 3 r c were 
used to construct relative count rate profiles, where the rel¬ 
ative count rate is defined for each crater by dividing each 
time series measurement by the mean count rate within r c 
of the crater centre. Each time series observation provides 
an estimate of the relative count rate at a given r/r c and 
these are stacked together for different crater subsets. 

4.1. Radial count rate variations 

The results in this subsection show how the neutron count 
rate varies as a function of sub-detector point distance to 
the crater centre. Stacked subsets of similar-sized highland 
craters are used, as are data for different neutron energy 
ranges. 

Figure 3 shows the mean stacked, count rate (c) profiles, 
with observations from each contributing crater normalised 
by the mean count rate measured from positions over that 
crater, c(< r c ). As the craters are very well sampled, had 
the count rates instead been normalised with respect to the 
count rate at r/r c ~ 3 the stacked profiles would only change 
by a radius-independent, vertical shift. The radii are nor¬ 
malised by the relevant crater radius. Points and errors on 
the mean profiles are shown only for the 40 < r c /km< 50 
case for clarity, but are of similar size for the other crater 
subsets. For both the thermal and epithermal profiles, a cen¬ 
tral ~ 1% enhancement in the neutron count rate is seen, 
with the count rate outside the crater being ~ 0.7% lower 
than the mean count rate measured over the crater. These 
features are about twice as pronounced as those in the corre¬ 
sponding fast neutron profiles, and are common to all crater 
samples with r c > 40 km. For smaller crater sizes, the fea¬ 
tures in the profiles decrease in amplitude. Given that the 
FWHM of the LP neutron detectors is approximately 45 km 
at an altitude of 30 km, one should expect that any features 
on smaller scales will be washed out. Also, if all craters 
have their radii either overestimated or underestimated by 
5% then the changes in the neutron count rate profiles are 
lower than 0.1%, so the results are robust to this level of 
systematic uncertainty in the crater radius determination. 

One reason why the fast neutron profile might be less vari¬ 
able than those in the lower energy ranges is that the tem¬ 
poral sampling is lower, with 32 s observations, rather than 
8 s. At an altitude of 30 km, LP travelled ~ 50 km during 
32 s. The effect of the resultant blurring of the profile can be 
estimated by degrading the sampling of the epithermal neu¬ 
tron time series. This is illustrated for the 40 < r c /km< 50 
craters by the dotted line in the central panel of Fig. 3. For 
craters that are at least this large, the different sampling 
has only a small effect. However, for smaller craters, where 
the distance travelled during an individual observation cor¬ 
responds to a larger r/r c , the suppression of features in the 
normalised count rate profile will be larger. Fig. 3 suggests 
that any features in the fast neutron profile for craters with 
r c < 40 km would be small anyway. 

Lunar Prospector neutron count rates are evidently af¬ 
fected, at the ~ ±1% level, by the detector position relative 
to craters on the surface, provided that the detector foot¬ 
print is small enough to allow it to ‘see’ the craters. At this 
point, it is worth briefly considering the count rate profiles 
produced by the LEND SETN and CSETN. The SETN is 
an “omni-directional” detector, albeit strapped to the side 
of a “collimator”, so in practice it has an energy-dependent 
anisotropic footprint. Given that it is viewing the surface 
from ~ 50 km altitude, the features seen by the LPNS should 
be stronger than those recovered by the SETN. This is ev¬ 
ident in Figure 4, which shows the count rate profiles mea¬ 
sured by the SETN for a range of different crater sizes. The 
features, while still significant, have been washed out, typi¬ 
cally decreasing the amplitude of the central peak by a factor 
of a few. 

Figure 5 shows the spacecraft background-corrected 
CSETN count rate profiles. After correction, the count rates 


are typically only ~ 2 per second, hence the large statisti¬ 
cal uncertainties. However, unlike the SETN profiles that 
show similar trends with crater size to the LPNS results, the 
CSETN profiles show no obvious trends or significant central 
bumps in the count rate. This is entirely consistent with the 
large CSETN detector footprint inferred by Teodoro et al. 
[2014], 

Having determined that the LPNS count rate varies sys¬ 
tematically with distance from crater centres, the question 
becomes what is responsible for this? One uninteresting 
possibility can be immediately discounted by recalculating 
the count rate profiles using the raw LPNS data. The fea¬ 
tures are similarly present in the raw data, implying that the 
data reduction process did not create them and they do re- 



r / r c 


Figure 4. Stacked, normalised radial LEND SETN neu¬ 
tron count rate profiles for observations at altitudes less 
than 60 km. The different colours correspond to different 
sized craters, as described in the caption for Figure 3. 



r/r c 

Figure 5. Stacked, normalised radial LEND CSETN 
neutron count rate profiles for observations at altitudes 
less than 60 km. The different colours correspond to dif¬ 
ferent sized craters, as described in the caption for Fig¬ 
ure 3. 
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Figure 6. Variation of the model count rate profile normalised to the mean count rate over the crater. 
Panel (a) shows the variation with crater depth, d, for a 45 km radius crater, with no infill or uplift, 
and the detector placed at an altitude of 30 km. These model craters are just spherical cap depressions. 
Panel (b) fixes d = 3.0 km and shows the effect of adding a flat central infill region out to a fraction Xi/r c 
of the crater radius. Fixing Xi = 0.6 r c and x e = 2.0 r c , panel (c) illustrates the variation of the count 
rate profile with crater uplift. Finally, panel (d) shows how the model profile varies with crater radius. 
At each different radius, the crater shapes are defined using the appropriate parameters in Table 1. The 
data points show the LPNS count rate profiles for the corresponding crater stacks. 


fleet something to do with the lunar surface. Compositional 
variation would not create almost identical features in the 
thermal and epithermal count rate profiles. Also, if mafic 
and magnesian central peaks [Cahill et al, 2009] were having 
an important impact on these profiles, then the thermal and 
fast neutron profiles should be anticorrelated, whereas they 
show qualitatively similar behaviour. The possibility that 
these profiles are the result of the geometrical configuration 
will be considered in the following section. 

5. A Simple Geometrical Model 

A model describing how topography affects the detected 
neutron flux is outlined in this section, in order to determine 
if this alone can explain the neutron count rate profile over 
craters. The predictions from this model are compared with 


the LPNS count rate profiles and the implications of this 
comparison are then discussed. 

5.1. The model 

The flux measured a distance r away from a patch of sur¬ 
face area dA emitting neutrons at a rate fo per unit area, 
with a detector an angle 9 away from the surface normal can 
be written as 

df(r, 9) = (2 + “ )/o cos 1+Q 9 <LA, (10) 

where a represents the effective beaming of the neutrons out 
of the surface resulting from the increase in neutron num¬ 
ber density with depth in the top ~ mean free path in the 
regolith [McKinney et al, 2006]. a « 0.5 provides a good 
match to the Monte Carlo neutron transport flat surface 
models of Lawrence et al. [2006] for the range of neutron 
energies detected by the LPNS. 
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The model for the total flux received by the detector in¬ 
volves integrating equation (10) over the lunar surface that 
is visible from the detector, assuming that the flux from a 
particular piece of surface is proportional to the incident 
cosmic ray flux. One complication is that when the surface 
includes concave craters, their walls can act to block parts of 
the crater interior from the detector’s view. If crater uplift 
is included then this effect can also extend to the exterior 
of the crater. The model accounts for this but does not, by 
default, allow neutrons emitted from the crater and imping¬ 
ing on the crater wall to be reemitted. This assumption will 
be considered further in section 5.3. 

In practice, the flux calculation can be more efficiently 
performed by partitioning the surface into different zones 
and using symmetries in the problem to avoid needing to do 
a two-dimensional numerical integration over the full visible 
surface. These zones are: 

1. the central infilled region of the crater, (x < Xi ), 

2. the constant radius of curvature crater walls, (*i < 
x < x c ), 

3. the outer uplifted slope, ( x c < x < * e ), 

4. the unperturbed surface beyond the outer uplifted 
slope, (x e < x). 

For the more interested reader, Appendix A contains spe¬ 
cific details of the calculations involved. 

5.2. Predictions of the model 

Figure 6 shows the neutron count rate profile from the 
model. The four panels show how the profile changes with 
(a) crater depth, (b) extent of the infill region, (c) amount 
of uplift, and (d) crater radius. In all cases the detector is 
placed at an altitude of 30 km and the effective beaming of 
neutrons is taken to be a = 0.5. The mean altitude for the 
LPNS observations being considered at altitudes less than 
45 km is ~ 31km. 

Panels (a) and (b) show that, for a 45 km radius crater, 
as the crater becomes deeper or the infill zone smaller, the 
central peak in count rate over the crater increases. This 
happens because these changes enhance the effect of the 
neutron beaming seen over the crater. If the beaming of 
neutrons is switched off, i.e. a = 0, and neutrons are al¬ 
lowed to be reemitted off the crater interior, then the model 
has c(r)/c(< r c ) = 1 for all r. 

Panel (c) of Fig. 6 shows how uplifting a 45 km crater and 
having a constant gradient outer slope that returns to the 
unperturbed surface at x e = 2.0 r c affects the neutron count 
rate profile. As the uplift increases, the outer uplifted slope 
focuses more neutrons onto the detector when it is over the 
crater exterior, leading to larger count rates at r/r c > 1. 

The variation of the model count rate profile with crater 
radius is shown in panel (d). Parameters for the model 
crater shapes are taken from the fits to the stacked LOLA 
topographical profiles, as listed in Table 1. For the 15 km 
radius craters, the central peak in the count rate profile 
occurs on scales too small for the 45 km FWHM of an omni¬ 
directional detector at an altitude of 30 km. Consequently, 
the profile looks almost flat. For larger craters, the central 
peak in neutron count rate becomes increasingly apparent 
as the instrumental FWHM corresponds to smaller r/r c . 
The simple geometrical model captures much of the central 
bump that is present in the data for the different crater 
sizes. However, more apparent is the failure to reproduce 
the LPNS results at r/r c > 2, where the model overpredicts 
the observed count rate by ~ 0.3%. 

The features of the comparison between model and LPNS 
neutron count rate profiles are common across the different 
crater sizes, in both the deep and shallow craters, and when 


the observations are split into high and low altitude sub¬ 
sets and the model is adjusted accordingly. In all cases, the 
model appears slightly to underestimate the count rate ob¬ 
served over the crater. Given that this provides the normal¬ 
isation for all count rates, a consequence is that the model 
overestimates the normalised count rate at large distances 
from the crater. 

One might wonder if the stacking process, used here to 
increase the statistical significance of the measured average 
neutron count rate profile features, might introduce system¬ 
atic effects. For instance, not all craters in a particular 
radius range have identical topographical profiles. If the 
neutron count rate profile features were especially sensitive 
to the deepest craters, which might have only a small im¬ 
pact on the average topographical profile, then the stacked 
count rate profile might not reflect changes in the average 
topography. However, as the features in the neutron count 
rate profiles are small and the model performs similarly well 
for subsets of craters selected by radius or depth, this pro¬ 
vides reassurance that such non-linearities are unimportant 
here. Consequently, it is evident empirically that the model 
based upon the average crater topography does encapsulate 
the important features that are responsible for giving rise 
to the stacked neutron count rate profile, and the stacking 
procedure is an appropriate way to perform this study. 

5.3. Improvements to the simple model 

The small difference between the model and LPNS neu¬ 
tron count rate profiles presumably arises due to an inap¬ 
propriate assumption in the simple geometrical model. In 
this section, the assumptions being made in the model will 
be varied to determine what is required in order to fit the 
data. 

There is no energy dependence in the model predictions, 
so the similarity between observed thermal and epithermal 
count rate profiles and how they differ from the fast neutron 



Figure 7. Variation of the best-fitting neutron leak¬ 
age beaming parameter, a, with neutron energy, E. The 
fit to the angular distribution of leakage neutrons deter¬ 
mined using Monte Carlo neutron transport simulations 
has the flux as a function of angle from the surface nor¬ 
mal, 9, proportional to cos 1+ “ 9. Statistical uncertainties 
are smaller than the circles representing the simulation 
results, and the vertical line shows the lower limit of 0.5 
MeV for LPNS fast neutrons. 
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profiles is immediately suggestive that there is an energy- 
dependent misassumption in the model. The assumed beam¬ 
ing factor a = 0.5 is relevant for thermal and epithermal 
neutrons, but for LPNS fast neutrons, with energies above 
0.5 MeV, the best-fitting a decreases, corresponding to less 
beaming. Figure 7 shows these results from fits to Monte 
Carlo neutron transport simulations. Furthermore, the sin¬ 
gle parameter power-law fit does not accurately model the 
angular distribution of emitted neutrons at the fastest ener¬ 
gies, with the actual distribution being less beamed normal 
to the surface. Accounting for the LPNS detector response 
and the incoming flux as a function of fast neutron energy 
suggests that an appropriate value for a for the model is 
probably in the range 0.4 — 0.45. This lessening of the 
beaming acts to suppress the size of the features in the fast 
neutron count rate profile, and goes roughly half way to ex¬ 
plaining the difference between the LPNS epithermal and 
fast neutron count rate profiles. 

Another possible effect that might reduce the fast neutron 
count rate profile features is that the emitted fast neutrons 
may have an angular distribution that retains some memory 
of the direction of the incoming cosmic ray that produced 
them. The model assumes that the emitted neutron flux de¬ 
pends only on the angle from the normal to the surface, and 
not the azimuthal angle. Within craters, if fast neutrons 
are more likely to be emitted in the forward direction with 
respect to the incoming cosmic rays, then this would prefer¬ 
entially aim them into the crater and thus slightly reduce the 
count rate measured over the crater. This is qualitatively 
consistent with the difference between the fast neutron count 
rate profiles and the thermal and epithermal ones. Given 
these difficulties in modelling the fast neutron emission, the 
fast neutron results will not be considered further. 

Figure 8 shows the results found for the epithermal neu¬ 
tron count rate profiles of 45 km-radius craters when various 
different model assumptions are made. The common theme 
in tweaking the model is the desire to increase the count 
rate observed over the crater relative to that observed out¬ 
side the crater. For instance, the blue curve results from 
allowing all neutrons emitted from within the crater and 
aimed at another part of the crater interior to be reemit¬ 
ted rather than absorbed. Details of this calculation are 
described in Appendix (A5). The difference between no ree¬ 
mission and complete reemission, which is presumably also 
unrealistic, amounts to less than 0.1% in the count rate, so 
is insufficient to make the model fit the data. 

Figure 7 suggests that a = 0.515, rather than 0.5 rep¬ 
resents the best description of the beaming of thermal and 
epithermal neutrons, but this change is too small to make 
a significant difference in the count rate profile. Increasing 
the amount that neutrons are beamed from the surface by 
changing a from 0.50 to 0.51 within the crater, while leav¬ 
ing a = 0.50 for the crater exterior, has a larger impact on 
the predicted count rate profile. This is shown by the green 
curve in Figure 8, but it still fails to fit the LPNS results. 
A similar result is found if the number of neutrons emit¬ 
ted per incident cosmic ray is increased by 0.5% within the 
crater only (red curve). While this approximately recovers 
the LPNS profile for r/r c > 2, it predicts a dip at r/r c ~ 1.2 
that is deeper than observed. 

No single parameter change that has been considered is 
able to recover the observed LPNS neutron count rate pro¬ 
files. However, a good fit can be found by including a combi¬ 
nation of 75% neutron reemission from the crater walls and 
a 0.35% enhancement in the neutron yield for the region 
out to 2r c , which includes the crater interior and most of 
the outer uplifted slope. This fits the LPNS count rate pro¬ 
files constructed from time series observations taken from 
altitudes below 45 km. The same model also fits the profiles 
observed in different altitude ranges, as shown in Figure 9. 
At larger altitudes, the detector footprint is larger, and this 



Figure 8. Variation of the model normalised neutron 
count rate profile with different changes to the default 
model for a 45 km-radius crater (black line) relative to 
the LPNS epithermal neutron result for the stack of 
40 < r c /km < 50 craters (points). The red curve shows 
the model resulting from increasing the emitted epither¬ 
mal neutron flux per input cosmic ray by 0.5% within 
the crater. Including reemission of all neutrons (rather 
than the default of none) from the crater interior sur¬ 
face changes the default model to that shown with a blue 
curve. The green line results from increasing the neutron 
beaming from the surface from 0.50 to 0.51 within the 
crater only. 

suppresses the amplitude of the features in the count rate 
profile, and this behaviour is accurately captured by the 
model. This success was not inevitable, but it was neces¬ 
sary if the tweaks to the model are to be interpreted as 
telling us something about the lunar surface. Had Figure 9 
included results from craters/basins with sizes comparable 
to the LPNS detector footprint at an altitude of 100 km, 
then the high-altitude data would show a central count rate 
peak and a drop outside the crater radius. 

Similar changes to the simple model are able to fit the 
thermal and epithermal neutron profiles for all crater size 
ranges. The required enhancement is 0.3 — 0.4% for the 
crater subsets with r c < 50 km, with this enhancement rang¬ 
ing out to 1.5 — 2r c from the crater centre. 

One possible explanation for the enhanced neutron emis¬ 
sion could be surface or near-surface roughness, of the sort 
seen by radar measurements out to twice the crater radius 
[Stickle et al, 2015]. While previous neutron transport sim¬ 
ulations for planetary surfaces have assumed emission from 
a flat surface [Lawrence et al., 2006], it has been shown that 
the neutron leakage flux can be enhanced for non-flat sur¬ 
faces [Druke and Schaal, 1991]. If such roughness leads to 
the increase in emitted neutron flux required to fit the ob¬ 
served LPNS count rate profiles, then the neutron count 
rate could actually be sensitive to the physical condition 
of the lunar surface, making it complementary to the radar 
and thermal infrared data sets [Bandfield et al., 2011; Ghent 
et al., 2015]. The impact of changing the regolith mass dis¬ 
tribution near the surface can be addressed directly using 
Monte Carlo neutron transport simulations that use realis¬ 
tic topographic models of a planetary surface, as has been 
done for other planetary bodies [Prettyman and Hendricks, 
2015]. 
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6. Implications for hydrogen in polar cold 
traps 

The beaming of neutrons increases the count rate mea¬ 
sured by the LPNS when it passes over craters. This is the 
opposite effect to that produced by placing hydrogen into 
permanently shaded regions (PSRs) within polar craters, 
which reduces the epithermal count rate. Not accounting 
for the varying topography will lead to underestimates of the 
water-equivalent hydrogen cold-trapped into polar PSRs. If 
the change in observed count rate due to topography is 
~ 1%, then one might wonder how this could possibly have 
a significant impact upon the inferred WEH. However, this 
small change to the observed count rate is evident over the 
entire crater area, whereas the PSR may only cover a tiny 
fraction of the crater area. The blurring caused by the re¬ 
sponse function of the LPNS can have the effect of levering 
a small effect acting over the large crater area into a large 
effect in the small PSR area. 

Rather than considering a general crater, it makes sense 
to focus on Cabeus, where LCROSS actually made a local 
estimate of the WEH weight percentage. Cabeus is also 
significantly deeper than the average crater with a similar 
radius. The azimuthally-averaged shape of Cabeus is best 
fitted with a profile defined by d = 6.05 km, x-,/x c = 0.7, 
Xe/x c = 2.3 and u = 0.85 km. Adopting 42 km as the ra¬ 
dius of the crater [Head et al, 2010], and choosing a central 
circular disc covering 275km 2 as the PSR [Teodoro et al, 
2010], a model where the PSR region emits 0.35 times as 
many epithermal neutrons per incoming cosmic ray as the 
rest of the surface gives rise to the black line in Figure 10. 
This reasonably fits the epithermal neutron data for Cabeus, 
shown with black filled circles. The thermal neutron data, 
in contrast, are well-fitted by a model where the surface in¬ 
terior to 2.2r c emits 0.97 times as many neutrons as the rest 
of the surface. 

That the thermal and epithermal count rate profiles differ 
is consistent with the suggestion that hydrogen in the PSR 
is responsible for the odd shape of the epithermal neutron 



r/r c 

Figure 9. Stacked, normalised radial LP epithermal neu¬ 
tron count rate profiles for 40 — 50 km radius craters as a 
function of LP altitude. Blue, black and red results corre¬ 
spond to detector altitudes of a < 33 km, [33,60] km and 
a > 90 km respectively. Curves show the corresponding 
best-fitting model results. 


profile, although the reason for the thermal neutron pro¬ 
file rising above 1 at r > 2r c is not clear. If one ascribes 
the lack of an epithermal central count rate bump entirely 
to hydrogen in the PSR, then the factor of 0.35 in neutron 
count rate can be converted using the formula supplied by 
Lawrence et al. [2006] into ~ 4.5 wt% WEH. This is a factor 
of four greater than was inferred by Teodoro et al. [2010], 
and consistent with the LCROSS results [ Colaprete et al., 
2010; Strycker et al., 2013]. As Cabeus does not have a sim¬ 
ple crater morphology and the possibility that additional 
compositional variation is being suggested by the thermal 
neutron profile, this value should be taken with a pinch of 
salt. However, it serves to illustrate how the effect of to¬ 
pography upon remotely-sensed neutron count rates could 
lead to a significant bias in the inferred wt% WEH. Given 
that many PSRs occupy a larger fraction of the area of the 
crater within which they reside, the case of Cabeus may be 
a more extreme example of how large an effect topography 
can play. 

7. Conclusions 

This study shows that there are some features in the 
neutron count rate profiles sensed from orbital detectors as 
they are flown over lunar craters located in highland regions. 
There is a central bump in the detected count rate, and the 
mean count rate over the stacked crater is up to 1% larger 
than it is outside. This factor is largest for thermal and ep¬ 
ithermal neutrons, but still detectable in the fast neutrons. 

A simple geometrical model has been developed. It pre¬ 
dicts qualitatively very similar behaviour to that observed 
from the LPNS thermal and epithermal data sets. The cen¬ 
tral peak results from the weak beaming of emitted neutrons 
normal to the surface [Lawrence et al., 2006], which is anal- 



Figure 10. LP epithermal (black filled circles) and ther¬ 
mal (red open circles) neutron count rate profiles for 
Cabeus crater. Curves show models for an r c = 42 km 
crater with the topography that best-fits that of Cabeus, 
and 75% reemission of neutrons hitting the crater inte¬ 
rior. The black curve additionally has a central “compo¬ 
sition” that emits only 0.35 times as many neutrons per 
cosmic ray as the rest of the surface. The red curve traces 
a model where the surface interior to r c = 2.2 emits 0.97 
times as many neutrons per cosmic ray as the rest of the 
surface. 
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ogous to solar limb darkening. This simple model under¬ 
estimates the mean count rate observed over the crater by 
~ 0.3%. 

To fit the observed stacked count rate profiles well re¬ 
quires a ~ 0.35% enhancement in the neutron emissivity of 
the regolith within ~ 2r c of the crater centre. It should be 
possible, using Monte Carlo neutron transport simulations, 
to determine if this can be achieved by a plausible amount 
of surface or near-surface roughness. 

The beaming of neutrons over polar craters hosting PSRs 
may mean that the concentration of hydrogen in the PSRs 
has been underestimated in previous work. For the particu¬ 
lar case of Cabeus, where a large crater contains a relatively 
small PSR, it was shown that ~ 4.5 wt% WEH within the 
PSR can reproduce the epithermal neutron count rate pro¬ 
file, assuming a simple azimuthally-symmetric topographical 
model for Cabeus. This is a factor of four times larger than 
previously inferred, and is consistent with the value mea¬ 
sured using LCROSS data. In polar craters where the PSR 
occupies a larger fraction of the crater, the impact of topog¬ 
raphy on the inferred wt% WEH will be less important. 
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Appendix A: Details of the model neutron 
count rate calculation 

In order to calculate the model neutron count rate in an 
orbiting omni-directional detector, it is easiest to split the 
surface into a few distinct regions: the unperturbed surface, 
the outer uplifted slope, the crater walls and the flat infill 
region in the crater centre. The simpler case with no uplift 
will be considered first. 

Al. Neutron count rate from the uncratered surface 

This part of the calculation is very similar to that de¬ 
scribed in Appendix B of Prettyman et al. [2006]. The flux 



Figure 11. Variables used in the calculation of the count 
rate from the unperturbed lunar surface. The origin is 
placed at the lunar centre, O, and the detector, D, is at 
an altitude a above the surface. 


D. 



Figure 12. Variables used in the calculation of the count 
rate from the unperturbed part of the surface when a 
crater is present. The dotted line represents the intersec¬ 
tion of the cone with half-opening angle /3 and the sphere 
with radius r m . 


of neutrons a distance r away from a patch of lunar surface 
of area dA, at an angle 9 to the surface normal, as shown in 
Fig. 11, will satisfy 

CO s 1+a 9 

d/(r, 9) oc -dA, (Al) 

where a represents the effective beaming of neutrons from 
the surface. Integrating over 2n steradians and defining the 


z A 

F l 


V' 



Figure 13. Variables used in the calculation of the sky 
visibility from a point P, at y = 0, within the spherical 
cap crater. Points F and Q represent the crater focus and 
the point on the crater rim that lies in the vertical plane 
through P making an angle <p with the x axis. # max (0) 
represents the maximum colatitude to which P sees cos¬ 
mic rays at an angle <j>. rj is the angle between the x axis 
and the point underneath Q in the x — y plane, such that 
tan?? = vq/xq. 
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flux through the surface as /o, leads to 


d/M) = 


(2 + a)f 0 


‘ 9 AA. 


2nr 2 

The flux detected from the whole surface is then 

(2 + a)fo 


(A2) 


/ = 


2?r 


d</> 


m 


— ) cos 1+Q <9sin/3d/3. 


(A3) 

As shown in Fig. 11, /? is the angle subtended at the lunar 
centre by the vectors to the detector and surface patch and 
/Imax = cos -1 [r m /(r m + a)] defines the lunar horizon for a 
detector at altitude a. Using the sine and cosine rules, 


. „ fr m +a\ . 
sin 9 *= ^ —-- J sm /? 


(A4) 


this /?: 


COS (\>rim 


COS 0 C en COS /? — COS 8 C 
sin Qcen sin/? 


(A10) 


A2. Cosmic ray occlusion within a crater 

The cosmic ray flux impinging upon a unit area of crater 
interior will be lower than that incident on the outside, con¬ 
vex surface. Under the assumption that the cosmic ray 
flux is isotropic, this is accounted for by replacing /o with 
/ofli/7r, where ST is the cos(incidence angle)-weighted solid 
angle of visible sky. At a point P within a spherical cap 
crater, this is given by 


fii = 





cos i sin #d#d <j>, 


(All) 


and 


r 2 = (r m + a) 2 + r 2 m - 2r m (r m + a) cos/3. (A5) 

Defining t = (r m + o)/r m yields 

/-p + «)A /*“ - (A6 , 

Jo (t 2 + 1 — 2t cos (?) 2 

which can be computed numerically to find the flux from 
the uncratered surface. 

The detector has been assumed to be omni-directional 
in the above calculation such that the detected count rate 
is merely proportional to the flux at the detector. The 
LPNS is, in fact, cylindrical and thus is not quite omni¬ 
directional. However, comparison of the inferred instrumen¬ 
tal point-spread function with that given by Maurice et al. 
[2004] shows them to be similar in shape to the extent that 
correcting for any differences has a negligible effect upon the 
results in this paper. 

When a crater is inserted into the surface, the integra¬ 
tion limits in equation (A3) need to be changed. If the 
crater centre lies at the spacecraft nadir, then the minimum 
/? is increased so that the integration starts at the edge of 
the crater. However, for a more general crater position it 
is necessary to find the range of azimuthal angle <j> that lies 
outside the crater as a function of /?. Figure 12 shows this 
more general configuration, where the crater centre subtends 
an angle 9 cen at the lunar centre. Without loss of general¬ 
ity, the detector and crater centre can both be placed in the 
x — z plane, where the axes have been chosen such that the 
y axis is into the page and the detector is placed on the z 
axis. The required </>(/?) can be found by determining the 
points where the ring of lunar surface at /? intersects with 
the plane containing the crater rim. Using the fact that the 
crater centre lies in the same plane as the crater rim, one 
can infer that the rim plane is given by 

x.n = r cen = r m cos 0 C , (A7) 


where the incidence angle, i, is the angle between the vector 
PF and the direction ( 9 , <j>) and 0 max is the maximum angle 
down from the z direction that lies above the crater rim, 
as shown in Figure 13. Q represents the point on the crater 
rim at this particular azimuthal angle, <f>, and the vector PQ 
makes an angle 0 max with the z direction. The incidence an¬ 
gle can be written in terms of the two angular coordinates 
as 


(r-c 


zp) cos 9 — xp cos 4> sin 9 

T curv 


(A12) 


where xp and zp are the x and 2 coordinates of point P. 

Redefining the origin of the coordinate system to be at 
the base of the crater, the position of Q is given by 


( x c cos rj \ 
x c sin r) , 

d J 


(A13) 


with rj being the angle between the x axis and the point 
beneath Q in the 2 = 0 plane. Point P has coordinates 


/ Xp = V curv sill \ 

x P = 0 J , 

\ Zp = T curv (1 - COSl/>) J 

(A14) 

where tj) is the angle between the -z direction and FP. The 
vertical plane containing P and Q has 

/ sin 4> \ 
h = — cos (j> , 

V 0 ) 

(A15) 

and satisfies 

x.n = xp sinc/i. 

(A16) 


Inserting Xq into this equation yields the following expres¬ 
sion for r/ as a function of (/>: 


where the unit normal to the plane is given by 



h = 



(AS) 


Using the fact that 


PQ.2 = COS0 max , 


(AIT) 


(A18) 


Noting the symmetry in the y direction and finding the so¬ 
lution when 

( sin /? cos <j> \ 

sin/?sin0 ) , (A9) 

cos/? J 

leads to the following expression for </>rim(/?), the angle that 
represents the fraction of n radians outside the crater for 


one can infer that 
COS $max 


■ 


yj 'x% + X 2 p — 2x c Xp COST/ + (d — Zp) 2 


(A19) 


For a choice of crater shape and distance from the crater 
axis, xp , equations (A17) and (A19) determine cos# ma x(</>), 
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which can then be used in conjunction with equations (All) 
and (A12) to determine the fraction of sky visible from this 
point within the crater. 

Extending this approach to the case where there is a flat 
infilled region in the crater centre is straightforward. In 
practice, a table of SI values as a function of xp is cre¬ 
ated once, and this is used, with interpolation, for the two- 
dimensional numerical integration to find the flux coming 
from within the crater. 

A3. Visibility of a surface patch from the detector 

Out to the lunar horizon the crater exterior is all visible 
to the detector in the case where there is no uplifted rim. 
However, there are parts of the crater interior that may not 
be visible to the detector. Consequently, it is necessary to 
see if the line-of-sight from the detector to the surface patch 
passes above or below the crater rim. 

Placing the origin of the coordinate system, O, at the 
lunar centre, and the detector at 

( sin 6d \ 

0 I , (A20) 

cos Od ) 

with the z axis going through the crater centre, a general 
point within the crater can be written as 


back down to the unperturbed surface is a new topograph¬ 
ical component that also suffers from these issues. In con¬ 
trast, the calculation of the flux from the crater itself is 
only slightly changed to account for the raising of the entire 
crater surface. 

A4.1. Cosmic ray occlusion 

Considering first the occlusion of cosmic rays from the 
outer uplifted slope and the unperturbed surface, the sym¬ 
metry is such that this is just a function of the distance from 
the crater centre. The outer uplifted slope is most conve¬ 
niently parametrised using an azimuthal angle, e and the 
fraction of the way down the slope from the rim to the un¬ 
perturbed surface, /. For a point P on the outer uplifted 
slope, the azimuthal variation of the maximum polar angle 
to which the sky can be seen, # max (<)>), will be set either by 
the unperturbed surface or the outer uplifted slope, depend¬ 
ing on which is hit first as the zenith angle increases. 

Choosing the z axis to pass through point P and the 
crater centre to lie in the x — z plane at x < 0 , the value of 
0max to the unperturbed surface is independent of <p. Simple 
trigonometry gives 


COS $max,l 



(A24) 


( r curv sin ip cos 7 \ 

I'curv sin ip sin 7 ) , (A21) 

r m cos 9 C — d + r curv ( 1 — cos ip) ) 

where 7 represents the angle around from the x axis to point 
P. The symmetry of the problem means that the contribu¬ 
tion to the flux coming from 0 < 7 < 7 r is the same as that 
from —iv < 7 < 0. Following a similar methodology to that 
adopted in Section A2, the normal to the plane containing 
O, P and D can be defined using n = x p x (x D — x p ). The 
point Q on the rim determining if the detector is above or 
below the crater rim as viewed from P can then be found as 
the solution to Xq.u = x P .n with an x coordinate between 
those of P and D. In this case, 


with rp being the distance of point P from the lunar cen¬ 
tre. For sufficiently extended outer slopes, it is possible for 
rp < r m , in which case cos# max ,i is set to 0 . 

It may be that the outer uplifted slope itself is the first 
piece of lunar surface to intersect the line-of-sight as the 
zenith angle is increased at a particular azimuthal angle. 
In this case, # max is set by the local slope at point P in 
the azimuthal direction, <p. For a small displacement on 
the uplifted slope having components da: and Ay, such that 
(ds ) 2 = (da :) 2 + (dj /) 2 and taxup =Ay/Ax, the maximum 
zenith angle to the outer uplifted slope can be found from 


cos 0 max ,2 ( 0 ) = sin 


tan 



(A25) 


( x c cos 7 
x c sin 77 
r m cos 9 C 


(A22) 


d z dz da; dz Ay 
ds dx ds dy ds ’ 


(A26) 


and the plane equation is used to determine 7 . For the de¬ 
tector to be able to see point P requires 


~ (—p ~ x P ) > . (x Q -x P ) 
-\x D -x p \ ~\ x Q -x p \' 


(A23) 


These equations, along with those from Sections A1 
and A2, allow the computation of the curves in the top two 
panels of Figure 6 . One and two dimensional numerical 
integrations are required to evaluate the flux from outside 
and within the crater respectively. For the flux from within 
the crater, it is also necessary to compute the distance to 
the detector and the angle between surface normal and the 
detector direction, but these are readily found from the vec¬ 
tors used to determine if the detector can see that point 
within the crater. The two dimensional integration to find 
the crater flux is simply done over an azimuthal angle rang¬ 
ing from 0 to 2tt and the angle from the focus to the crater 
centre, ip, running from 0 to t/Wx = tan _ 1 [a; c /(rcurv — d)]. 

A4. Uplifted crater rim 


with dz/dy = 0 , ds/da: = l/cos<)> and dz/dx being the 
gradient g from equation ( 6 ) rotated through 0 cen into the 
coordinate system with P on the z axis. This leads to 

dz _ sin 9 cen + g cos 9 cen (A27) 

dx cos 9 cen — g sin 9 cen ' 

The value of cos # ma x(0) is taken as the larger of cosHmztx.i 
and cos 0 max , 2 ( 0 ), and the cosmic ray occlusion factor, Qi/n, 
at a given fraction of the way down the outer uplifted slope 
is calculated using equation (All). 

The cosmic ray occlusion for points on the unperturbed 
surface, like that on the outer uplifted slope, is just a func¬ 
tion of distance to the crater centre. It is convenient to 
place the patch of unperturbed surface under consideration, 
P, on the a axis and rotate the crater centre through an 
angle —9 cen about the y axis (moving the crater in the —x 
direction). Points an angle t around from the x axis on the 
crater rim, x R , or the outer edge of the outer uplifted slope, 
x E , can then be described via 


Including an uplifted crater rim complicates the calcu¬ 
lation considerably, because parts of the previously unper¬ 
turbed surface may now undergo some cosmic ray shadowing 
and may also no longer be visible from the detector. Sim¬ 
ilarly, the outer uplifted slope going from the crater rim 


( —(rm cos 9 C + u) sin 9 cen + rm sin 9 C cos e cos 9 cen \ 
r m sin 9 C sin e j 

(r m cos 9 C + u) cos 9 cen + r m sin 9 C cos e sin 9 cen ) 

(A28) 
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and 

( — r m cos 9 e sin 9 cen + r m sin 6 e cos e cos 9 cen \ 
r m sin 9 e sin t j 

r m cos 9 e cos Seen + r m sin 9 e cos e sin 0 cen / 

(A29) 

respectively. For a given azimuthal angle <f> around from the 
x axis as seen from point P on the z axis, it is possible to 
find any points on the uplifted outer slope that lie in the 
plane 


( sin <j> \ 

- cos 0 1=0. (A30) 


If there are no solutions for / in the range [0,1], then the 
plane at fixed <j> does not intersect the uplifted region and 
cos 0 max (<?i) = 0. If solutions exist, then equation (A30) pro¬ 
vides a constraint on /(e) for points on the outer uplifted 
slope that lie in the plane an azimuthal angle <j> around from 
the x axis as viewed from point P on the unperturbed sur¬ 
face. The largest zenith angle from which cosmic rays arrive 
at point P, 0max(0), is found using a numerical minimisation 
algorithm applied to the set of points on the slope. A root¬ 
finding algorithm is employed to determine / at any given 
e as part of this process. Given cos 0 max (</>) as a function 
of distance from the crater centre, the cosmic ray occlusion 
factors can be found using equation (All). 

A4.2. Visibility of surface patch from detector 

Points on either the outer uplifted slope or the unper¬ 
turbed surface may not be visible from the detector as a 
result of the uplifted region surrounding the crater. 

For a point P on the outer uplifted slope to be visible 
from the detector, D, the line of sight must not be blocked 
by either the outer uplifted slope or the unperturbed sur¬ 
face. If the dot product of the surface normal at P and 
the surface-to-detector vector, A = x D — x p , is positive, 
then P is not blocked by the outer uplifted slope. The 
unperturbed surface will block the line-of-sight if the line 
connecting P to the detector passes within r m of the lunar 
centre. Defining the fractional distance along this line as v, 
such that x = x p + vA, this happens if 0 < t> m i n < 1 and 
|*(vmin)| < r m , where u m i n represents the v for which this 
line passes nearest to the lunar centre. With the coordinate 
system origin at the lunar centre, this leads to 


v 


min — 


x P . A 

w 


(A31) 


and 


|*E(t , min) | 


2 


r 


2 

m 


\x P \ 2 - (x P .Af 
r 2 


(A32) 


where A is a unit vector in the direction of A. These equa¬ 
tions allow a quick determination of whether or not the un¬ 
perturbed surface blocks the detector’s view of a part of the 
outer uplifted slope. 

To determine the visibility of the unperturbed surface 
from the detector, consider placing the detector on the a 
axis at ( 0 , 0 , r,„ + a) and the crater centre in the x — z plane 
at x < 0. If the far point of the edge of the crater outer 
uplifted slope is visible above the far point of the rim (e = 7 r 
in equations (A29) and (A28) respectively), then the entire 
unperturbed surface is visible from the detector. If this is 
not the case, then the plane containing the lunar centre, de¬ 
tector and point P can be found. The line of intersection of 
this plane with the uplifted slope and the minimum zenith 
angle from P to points on this line follow, and the visibility is 
determined by comparison with the zenith angle from point 
P to the detector. This is a very similar methodology to 
that described to determine the cosmic ray occlusion factor 
for the unperturbed surface. 

A5. Neutron flux impinging upon the crater walls 

In the preceding sections of this appendix, the assump¬ 
tion has been made that any neutrons emitted from within 
the crater and aimed at the crater walls are absorbed on 
contact with the regolith and do not contribute to the neu¬ 
tron flux emerging from the crater. This is a simplification, 
because some of these neutrons will be re-emitted before be¬ 
ing absorbed. The more energetic neutrons may even lead 
to nuclear reactions that create more than one lower en¬ 
ergy neutron that escapes from the crater, in which case the 
crater would be producing a thermal neutron flux that was 
an amplified version of the incident fast neutron flux. 

While quantifying the impact of this process requires 
Monte Carlo neutron transport simulations, it is possible 
to use the simple model to determine how much of the emit¬ 
ted crater flux impinges on the crater surface as a function 
of position within the crater. Following the methodology of 
the previous sections, when the crater flux at the detector 
was determined, it is possible to place the ‘detector’ on the 
crater surface and calculate the flux from the crater that 
is aimed into the crater surface. The only additional fac¬ 
tor to consider is to include the fact that the normal to the 
‘detecting’ surface is at different angles to the lines-of-sight 
to the various other bits of crater surface. Multiplying the 
detected flux by the cosine of the incidence angle and in¬ 
tegrating over the entire crater surface leads to the results 
shown in Fig. 8 for the case where all neutrons are assumed 
to be reemitted from the surface. 



